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Abstract 

We show that, for Majorana neutrinos propagating in a constant magnetic 
field, the flavour survival probabilities for left-handed neutrinos is the same 
as for right-handed neutrinos, i.e., P^\vaL T^ah) = P^^{i^aR '-'aR), 
where a = e, //, r, whereas in the Dirac case the corresponding probabilities 
P^i^aL T^ah) and P^{h'aR ^aR) are in general different. This might 
lead to a novel way to search for the nature of neutrinos. We also discuss how 
this relation for Majorana neutrinos gets modified when the magnetic field is 
not constant. However, if matter effects become important the relation does 
not hold anymore. 
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I. INTRODUCTION 



In the past few years neutrino physics has generated tremendous interest and has been 
looked upon as a window to physics beyond the standard model. Recent results from Super- 
Kamiokande |]T| have boosted this search for new physics. Both the solar and atmospheric 
neutrino anomalies agree well with the neutrino oscillation hypothesis (for fits to the 
solar neutrino data see Ref. |^), which requires neutrinos to be massive. For recent reviews 
see, e.g., Ref. [0. However, it should be kept in mind that neutrino flavour oscillations is 
not the only possibility to describe the data. It may be that massive neutrinos have non- 
zero magnetic moments (MM) and electric dipole moments (EDM), which could play an 
important role in the solar neutrino puzzle if large magnetic fields are present in the interior 
of the sun ^ . In this framework, the attractive scenario of resonant spin - flavour transitions 
has received attention and good fits to the solar neutrino data have been obtained (for 
recent fits see Refs. ^ and for reviews see Ref. 0]). In solutions of the solar neutrino puzzle 
with non-zero MMs and EDMs, it is assumed that the solar neutrinos interact with the 
magnetic field in the sun to produce right-handed neutrinos due to a helicity flip. Right- 
handed neutrinos are sterile in the case of Dirac neutrinos and behave like antineutrinos in 
the Majorana case. In the latter case, a flavour transition simultaneous with the helicity 
flip produces the suppression in solar neutrino detection with elastic neutrino - electron 
scattering. 

So far, it is not known if neutrinos possess MMs and/or EDMs. The most stringent 
laboratory bounds come from elastic neutrino - electron scattering: for reactor neutrinos 



the limit is 1.8 x 10 ^"/is whereas for the limit is 7.4 x 10 ^V-B [0) where hb is the 



Bohr magneton. More stringent limits are obtained with astrophysical considerations [ITT 
For a collection of references on limits on neutrino MMs see also Ref. [jl2| 



In this paper we focus on neutrino survival probabilities. Apart from the three active 
neutrino flavours (z/^ with a = e, h,t), we allow an arbitrary number of additional neutrinos 
of the sterile type [ug) ||13|. In vacuum, the survival probabilities of left-handed and right- 



handed (anti) neutrinos are equal: 

where the superscripts D and M refer to Dirac and Majorana neutrinos, respectively. These 



equalities are not valid if matter effects [|14[ become important, but matter effects in neu- 
trino oscillations do not distinguish between the Dirac and Majorana nature [|T^]. Here 
we concentrate on the situation that matter effects are negligible but the MM and EDM 
interaction of neutrinos with magnetic flelds becomes important. Thus our discussion does 
not apply to solar neutrinos. We will show that for constant magnetic flelds the equality 
(H) remains valid for Majorana neutrinos whereas in general it gets lost for Dirac neutrinos. 
This is to be contrasted with the matter effects. We will also consider a situation where the 
magnetic fleld depends on x where this equality for Majorana neutrinos persists. Though our 
discussion will be rather formal we hope that it can shed some light on a possible distinction 
between Dirac and Majorana neutrinos. Up to now such efforts are mainly concentrated on 
neutrinoless double beta decay p6[] . 

Our paper is organized as follows. In Section II, we discuss the general formalism for 
an oscillating neutrino (both, left and right-handed neutrinos) propagating in an external 
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electromagnetic field, in the absence of matter effects, and prove the equality ([l|) for Majo- 
rana neutrinos in a constant magnetic field. In Section III we study some conditions where 
equality (|ID also holds for Dirac neutrinos, but we show that it is violated in the general 
Dirac case. In Section IV we consider non-constant magnetic fields and in Section V we 
summarize our results. 



II. GENERAL FORMALISM 

In the following we will work with n neutrino flavours or types. We will take into account 
neutrino mixing and electromagnetic interactions through MMs and EDMs but we will not 
consider matter effects. 



A. DIRAC CASE 

Let us assume that the neutrino MMs, EDMs and transition MMs and EDMs are given 
for the neutrino mass eigenfields Vj. Then the MM and EDM interaction of the Dirac 
neutrinos is expressed by the Hamiltonian density 

= l^if^ + id-i^)ao.pF''^y with = = d (2) 

being the magnetic moment and electric dipole moment matrices, respectively. F""^ is the 
antisymmetric electromagnetic field tensor. Whereas Vj denotes the fields in the mass basis, 
we denote the chiral fields in the flavour basis where the charged lepton matrix is diagonal 
by ur. The mass matrix M in the neutrino mass term 



is bidiagonalized by 



CZ = ^RMiyL + h.c. (3) 
U^MUl = M (4) 



such that 

^aL = ^Laj VjL and V^R = ^ URaj l^jR , (5) 

j j 

where the indices a denote the neutrino flavours or types. The diagonalizing matrix Ul is 
the usual neutrino mixing matrix. 

For massive Dirac neutrinos with mixing, the Hamiltonian in the flavour basis describing 
the neutrino interacting with a magnetic field is given as [|^,|T7|,|r^ 

^ \-B_UR{^^-^d)Ul ^UrM^U^ j' 

where the upper half of the matrix corresponds to negative helicity while the lower 
half corresponds to positive helicity. Assuming that the neutrino is propagating along the z 
direction, the magnetic fields B± are defined as 
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Bx i iBy 



Be^'f^ with B 



Bl + Bl 



(7) 



In the approximation we are working with, the longitudinal magnetic field is negligible. The 
neutrino energy is denoted by E^. Note that for the right-handed Dirac neutrino states there 
is no preferred fiavour basis because these states do not couple to the charged leptons. The 
Hamiltonian matrix H? for Dirac antineutrinos is obtained by making the replacements 



U 



L,R 



d 



-d* 



(8) 



in the matrix (^. The superscript * on the MM and EDM matrices indicates complex 
conjugation of all elements of the matrices. The upper half of HP maps onto positive and 
the lower half onto negative helicity. 

The Hamiltonian (^) is re-expressed in the mass basis as H^,, where 



(9) 



and the corresponding Hamiltonian for antineutrinos is given by 



B4fi 



B+{fi* + id* 
id*) T^M^ 



(10) 



2E^ 



In order to obtain the survival probabilities of neutrinos (which have negative helicity) 
and antineutrinos (which have positive helicity), one needs to diagonalize the matrices H^, 
and Hp, respectively. Observing that 



J^H-J = H* with J 



1 
-1 



(11) 



where 1 is the n x n unit matrix, we find that H^ and ifp have the same eigenvalues 
El, ... , E2n- Furthermore, there are unitary matrices and Wp such that 

WlH^W, = WlHpWp = diag(Ei, ^2, E2n) , (12) 

which, according to Eq.(|n|), are related by 

w, = jw: . (13) 

-1) and Dirac antineutrinos 



The survival probabilities for Dirac neutrinos (helicity 
(helicity = -|-1) are then given as 



aL 



2n 

Hp. 



Oi] 



^-iEjL 



(14) 



and 



aR 



l^aR 



2n 

EP 



(15) 
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respectively. In the above probability expressions, L is the distance between neutrino source 
and detector. From the 2n x 2n matrices If^ and U" diagonalizing Hj^ and HP, respectively, 
we only need the first n lines labelled by the neutrino flavours or types, which can be 
expressed by f/^ and W = W^: 



l<J<n: = ELi ULckWk, , = ELi Ut^,W:^,^ , 

n + l<j<2n: U^^ = ELi ^La.W^.n+, , U'., = ELi t^£«.^nVn+, • ^ ^ 



B. MAJORANA CASE 

For Majorana neutrinos we start with the same electromagnetic Hamiltonian density 
(0) as in the Dirac case, except that the factor 1/2 is replaced by 1/4 to account for the 
charge conjugation property {I'jY = Vj of the Majorana fields. Then the Hamiltonian Hm, 
corresponding to , is given by 

Hm = 



^-UlM^uI -B+Ul{^ + id)Ul 



2E^ 



B^Ul{n-id)Ul 



(17) 



Note that Ur in is replaced by in the Majorana case. Furthermore, we have to 

keep in mind that the matrices /i and d are antisymmetric, i.e.. 



jj, , d^ = d* = —d 



(18) 



in Hm- This equation formulates the well known properties of Majorana neutrinos that only 
transition moments are non-zero and that /i and d are purely imaginary. The latter property 
allows to check immediately that 



J^HmJ = HZ 



M 



(19) 



holds (compare with Eq.(p])). 

Denoting the diagonalizing 2nx2n unitary matrix of Hm by Um, the survival probabilities 
of Majorana neutrinos with negative and positive helicities are given by 



In 



(20) 



and 



In 

E\Tj a -iEjL 



(21) 



respectively. The subscript — of Um in Eq. (pOD indicates that the index a refers to the first 
n lines of this matrix (negative helicity), whereas in Eq.(|21|) the subscript + indicates the 
use of the last n lines (positive helicity). 

The relation (|TP|) allows to be more specific with respect to the matrix Um- 
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Lemma: The matrix Um which diagonahzes Hm can be chosen to be of the form 

Um = 




Furthermore, diagonahzing Hm with this Um we obtain 



Ei = E, 



n+j 



for j 



n. 



Proof of the lemma: Suppose we have an eigenvector ip of Hm-, i-e., 

Hm^ = . 



Then with Eq.(|T9D we find that 



HmJ^P* = J {J^HmJ) ^* = JHljtjj* = EJtfj* and Jijj* ± ^ ^ 



(22) 



(23) 



This observation allows to construct an orthonormal basis of eigenvectors of Hm in the 
following way. Starting with an eigenvector ipi normalized to unit length and with eigenvalue 
El, then Jipl is orthogonal to ipi and has the same eigenvalue. Next we can find a normalized 
eigenvector ip2 with eigenvalue E2 such that ip2 -L {i^i, Ji^i}- Then it is easy to show that 
Jip2 is orthogonal to all three previously constructed eigenvectors. We continue by finding 
ips with eigenvalue E3 orthogonal to all four previously constructed eigenvectors and so 
on. After n steps, two orthonormal systems „ and {J4'j}j=i,...,n with the same 

eigenvalues Ei, . . . , E^ are found which together form an orthonormal basis of eigenvectors 
of Hm- Thus Um is given by 



UM = {i^i---i^nJri---jrn) 



(24) 



which is of the form announced in the lemma. □ 

An immediate consequence of the lemma and of Eqs. (pOD and (121]) is the following theo- 
rem. 

Theorem 1: Without matter effects and in a constant magnetic field the survival proba- 
bilities for left-handed and right-handed Majorana neutrinos are equal, i.e.. 



III. ASYMMETRIES OF SURVIVAL PROBABILITIES 

Let us define an asymmetry 

Af = P^'iv^L -> l^aL) - P^'i^aR ^aR) (25) 
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for every Dirac neutrino flavour or type a. Note that we have just discussed in Theorem 1 
that the corresponding asymmetry for Majorana neutrinos, 



A^^ = P'' {UaL ^ lyaL) - P"' {y^R - ^^ai?) , (26) 

vanishes with the assumptions in this theorem. If we set n = d = and assume vanishing 
matter effects, then we have vacuum oscillations and both asymmetries (P^D and (EBI) are 



zero. Clearly, the asymmetry (p5D is not a CP asymmetry due to a KM type of phase [jT9 
If it is non-zero it is because the presence of a magnetic fleld represents a CP-violating 
situation. 

As a passing remark, matter effects would induce such an asymmetry for the same reason: 
background matter is not CP-invariant or, in other words, neutrinos and antineutrinos inter- 
act differently with matter. The matter potentials enter in the diagonal of the Hamiltonian 
matrices as follows: 

(Vl, 0) for Dirac neutrinos, 
{—Vl, 0) for Dirac antineutrinos, (27) 
(Vl, —Vl) for Majorana neutrinos, 

with the diagonal matrix Vl = V^Gp diag(Na), where, for ordinary matter, N^. = ne — nn/2, 
Nfj^^r = —nn/2, Ns = and Ue and Un are the electron and neutron densities, respectively. 

However, there is a fundamental difference between matter effects and the effects of MMs, 
EDMs and a magnetic fleld: For Bj- = 0, but matter effects becoming important, one has 
= A^^ 7^ in general, whereas with Vl = but Bj- ^ and constant magnetic fleld one 
has A^ 7^ in general, as we will see, but A^^ = 0. 

In view of Theorem 1 and in order to elaborate under which conditions Majorana and 



Dirac neutrinos behave differently, it is necessary to study the asymmetry (|25D in more 
detail. The following theorem describes under which sufficient conditions the asymmetry is 
still zero. 

Theorem 2: For the case of Dirac neutrinos, without matter effects and with a constant 
magnetic fleld, if the electric and magnetic moment matrices obey the proportionality fi = cd 
(or d = cn), where c is a real number, then the survival probabilities of a neutrino with 
flavour (type) a is equal to the survival probability of the corresponding antineutrino, i.e., 
the asymmetry ([25| ) is zero. 

Corollary: If(i = OoryU = Oorn = l the asymmetry (|25|) is zero. 
Proof of the theorem: We deflne (see Eq.(|^)) 

e'^{l + ic)=z=\z\e'^. (28) 

We can write Hu @ as 

The phase ( can be removed by the unitary transformation 
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1 

le'^ 



1 

le~'^ 



-B\z\fi 





' Hi 


-H2\ 


M2 ] \ 


.-H2 


Hi 1 



(30) 



The individual block matrices Hi and H2 are independent Hermitian matrices. 
We consider the eigenvector equations 



(31) 



{Hi - H2)Xj = EjXj , 
(Hi + H2)Yj = E^^^Y.j , 

where j = 1, . . . ,n. The sets {-^j}j=i,...,n and {l^}j=i,...,n form orthonormal bases of eigen- 
vectors of the Hermitian matrices H1—H2 and H1+H2, respectively. Therefore, the structure 
of Wi, is given by 



Xi 

V2 [ e-'^-Xi 



Xn 

e~'^Xr. 



Yi ■■■ 

-e-^Yi ■■■ 



Y 



C 



D 



-e-^Yr,. 



(32) 



Furthermore, using the relations ([T6|) and the expressions (|T^ and ([T^) for the Dirac neutrino 
survival probabilities. Theorem 2 follows. □ 

Now we want to show that Theorem 2 describes an exceptional situation and that indeed 
in general the asymmetry is different from zero. It is sufficient to see this in the case of 
vanishing neutrino masses in H^^ (^) and H? (H). Defining a matrix [jl8| 



\ = ji — id . 



(33) 



we notice that this is a completely general matrix which can be bidiagonalized with unitary 
matrices R and S: 



A = R\S^ with = (xi, . . . , x„) , R= {yi, ...,yn) 



(34) 



where A is diagonal and positive and {xj}j=i^,,,^„ and {yj}j=i^...,n are orthonormal bases. 
Dropping the neutrino masses we get the Hamiltonian matrix 



i^.^ = -B 



e"^SXR^ \ 
e-'^RXS^ j 



(35) 



which has the following eigenvectors: 

1 



Xj 



v/2 V e-'^yj 



V2 



-e 'f^yj 



with eigenvalue —BXj , 
with eigenvalue BXj . 



(36) 



These eigenvectors of the Hamiltonian matrix (|35D form the matrix (p!2D, and thus with 
Eqs.([l3), dH) and (|16D we obtain 

2 



aL 



aR 



E 



X 



a] I 



cos{BXjL) 



\yaj\^ cos{BXjL) 



(37) 
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From these expressions it is obvious that in general the probabihties P^iyaL ^ah) and 
P'^i.^aR ~^ ^aii) are different, because the orthonormal bases {a;j}j=i^...^„ and {yj}j^i^ n are 
independent of each other. 

To elaborate this in more detail we consider the case of two neutrino flavours (n = 2). 
Since we neglect here the neutrino masses all phases in A ( |33D can be removed and the 
matrices S and R (plf) are characterized by the angles 9 and 9' , respectively: 



/ cos 6* - sin 6* \ , , / cos 6*' - sin 6*' \ , , , . 

\ sm 9 cos 9 j \ sm cos 9 1 ^ ' 

Inserting Xj and yj into the survival probabilities (0), it is evident that P'^iyai ^ai) 7^ 
P^ipaR ^aR) hoMs as loug as Ai 7^ A2 and cos^ 9 7^ cos^ 9' . 



IV. MAJORANA NEUTRINOS AND z-DEPENDENT MAGNETIC FIELDS 



Even if we assume that = but B^, By depend on the survival probabilities for 
left-handed and right-handed Majorana neutrinos will in general be different, i.e.. Theorem 
1 will not hold anymore. This provides some obstacle for an application of the results of this 
paper with the aim to find a way to distinguish between the Dirac and Majorana nature of 
neutrinos. We briefly discuss the general formalism for the case of z-dependent Hm- One 
needs to find 2n linearly independent solutions of the differential equation 



.dip{z) 
dz 



Hm{z)>^{z) 



(39) 



The survival probabilities depend now also on the locations zq of the neutrino source and zi 
of the neutrino detection. Assuming to know a complete orthonormal set {v5j(2;)}j=i,...,2n of 
solutions of Eq. (|39D , we can formulate the transition and survival probabilities as 



P''\yaR{zo) ^ ypR{z,)) 



2n 



2n 



II'^+aj(^0)'^+/3i(2;i' 



(40) 



and similarly those with LR and RL transitions. We use the indices —a and +a to indicate 
the n upper and n lower components of <y9, respectively. 

If the matter potential can be neglected (V^ = 0), the relation (^) for the Majo- 
rana Hamiltonian matrix holds, irrespective if Hm depends on z or not. It is then easy 
to check that, given a complete orthonormal system of solutions {v?i(-2)}j=i,...,2n., then 
{ JV9*(— z)}j=i,...,2n is a complete orthonormal set of solutions of the differential equation 
(p9|) with Hm{z) replaced by Hm{—z). This suggests that with a Hamiltonian matrix Hm 
symmetric in z one can arrive at a situation where Theorem 1 is still valid. Indeed, using 
Eq.(^Dj) and replacing (pj{z) by Jlp*(—z), one arrives at the following relations for the Ma- 
jorana transition and survival probabilities. 
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Theorem 3: 

Vl = 0, HMiz) = Hm{-z) ^ P'"\vc.l{zo) ~^ mi^i)) = P^\^pr{-zi) ^ Vo.r{-z^)) , 

(41) 

and, therefore, with zi = —zq, we have 

P^'{iyaL{-Zi) ^ U^Lizi)) = P''{iyaR{-Zl) U^Zi)) . (42) 

Therefore, Theorem 1 holds also for a ^-dependent magnetic field, provided it is symmetric 
between the neutrino source and detection points. 

Note that from Eg . (1401) by exchanging complex conjugation within the absolute values, 
for any Vl and z-dependent magnetic field, we obtain the relations 

P^'{yaL{z^) ^ = P^'imi^l) ^ ^aL(^o)) , 

P'^'ii^aRizo) ^ Ufsniz,)) = P''Mz,) u^nizo)) , (43) 

which express CPT invariance. However, they do not allow to relate the survival probabilities 
for different helicities and are, therefore, not useful in our context. 



V. SUMMARY 

In this paper we have studied the survival probabilities of neutrinos and antineutrinos 
which possess magnetic moments and electric dipole moments and propagate in magnetic 
fields. In particular, given a neutrino flavour (type) a, we have studied the difference 

Af*'^^ (|25| , |26D between the survival probabilities for UaL and 'z^'aj?. The bar on UaR refers 
to Dirac neutrinos, without bar a Majorana neutrino is understood. It is well known that 
matter effects lead to non-zero asymmetries A^'*^, but without electromagnetic neutrino 
interactions one has A^ = A*^ expressing the fact that with neutrino oscillations one cannot 
distinguish between Dirac and Majorana neutrinos. 

Here we have considered the opposite situation: we have assumed that matter effects are 
negligible but neutrinos have MMs and EDMs and propagate in magnetic fields. Assuming 
a constant magnetic field, we have shown that in this case we have A*^ = (Theorem 1), 
but A^ 7^ in general. However, if the MM and EDM matrices are proportional to each 
other in the case of Dirac neutrinos, then the asymmetry of survival probabilities is zero as 
well, i.e., A^ = (Theorem 2). What happens if the magnetic field is not constant along 
the neutrino path between source and detector? In general, this will result in A^^ 7^ for 
Majorana neutrinos, but if the transverse magnetic field is symmetric with respect to the 
center Zc of the line connecting source and detector {B^^y^Zc + z) = B^^y^z^ — z)) one still 
has A^^ = (Theorem 3). 

The observations summarized here indicate a fundamental difference in the behaviour 
of Dirac and Majorana neutrinos with respect to magnetic fields, which is a consequence of 
the fact that the Hermitian MM and EDM matrices /i and d have to be antisymmetric for 
Majorana neutrinos. Further research is necessary to see if the results of this paper lead to 
realistic possibilities to distinguish between the Dirac and Majorana nature of neutrinos. 
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